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Abstract 

We revisit the two vortex merger problem (both symmetric and asymmetric) for the Navier-Stokes 
equations using the core growth model for vorticity evolution coupled with the passive particle field and 
an appropriately chosen time-dependent rotating reference frame. Using the combined tools of analyzing 
the topology of the streamline patterns along with careful tracking of passive fields, we highlight the 
key features of the stages of evolution of vortex merger, pinpointing deficiencies in the low-dimensional 
model with respect to similar experimental/numerical studies. The model, however, reveals a far richer 
and delicate sequence of topological bifurcations than has previously been discussed in the literature for 
this problem, and at the same time points the way towards a method of improving the model. 



1 Introduction 



1.1 General background 

Consider an incompressible fluid in an unbounded two-dimensional domain R 2 . The fluid motion is governed 
by Navier-Stokes equations, written in terms of the vorticity cj(x, t) as follows 



duo 



- — u • Vcj + vAuj , 



(i) 



where v is the kinematic viscosity taken to be constant, u(x, t) is the fluid velocity function of time t and the 
position vector x. Both x and u can be expressed in inertial frame {e^ = i 5 2,3 where (ei,e2) span the plane 
of motion. For simplicity, we denote the nonzero components x = (x, y) and u = (u Xl u y ). The vorticity 
vector uj = V x u which is always perpendicular to the plane of motion is therefore denoted by its component 
uj in the e3-direction. The velocity u and vorticity uj are related via the 2D Biot-Savart law 
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where x is an integration variable and x x = (—y,x). The solution of the system of equations ([!]) and Q 
depends on the initial condition cj(x, 0). One solution of particular interest in this work is the Lamb- 
Oseen solution corresponding to a point vortex initial condition cj(x, 0) = r<5(x), centered at the origin 
with strength T. Traditionally, the problem can be expressed compactly in complex notation with position 
variable z = x + iy. The Lamb-Oseen vorticity and velocity fields are given by 
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where the notation () = d()/dt refers to the time derivative and the notation z* refers to the complex 
conjugate of z. The vorticity field is initially concentrated at the origin and then diffuses axisymmetrically 
as a Gaussian distribution. The spreading of the vorticity can be quantified by the vortex core defined as 
a = \[Avt \/At. 

For more complicated initial data, exact solutions of ([!]) and Q with "growing cores" are not an- 
alytically available in general. We are particularly interested in the viscous evolution of a class of initial 
vorticity fields cj(z,0) = J2a=i ^0^(2 — z a ) consisting of the superposition of N point vortices. Gallay & 
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Fi gure 1: Schematics of (a) symmetric and (b) asymmetric co-rotating vortex pairs. 



Wayne proved recently that the Lamb-Oseen solution is an asymptotically stable attracting solution for all 
(integrable) initial vorticity fields, see pQ, of which the N point vortices is a special case. Our goal in this 
paper is to elucidate the dynamics that takes place as a given system evolves to the Lamb-Oseen solution. 
We use a familiar "core-growth" model (see [2]) on the two- vortex merger problem in an appropriate rotating 
reference frame with the knowledge that this model is not exact, but is nonetheless extremely useful when 
carefully implemented on a problem where the intermediate evolution has been carefully described both ex- 
perimentally and numerically in the literature. The comparisons will allow us to pinpoint the shortcomings of 
the model in this context and build "higher fidelity" , but still low-dimensional physically based corrections. 

We examine the viscous evolution associated with two initial configurations: a symmetric and an 
asymmetric co-rotating vortex pair. In other words, we consider two initial conditions: (i) two point vortices 
with equal circulation Tl = Tr = Y located at zl(0) = —bo/ 2 and z#(0) = bo/2 and (ii) two point vortices 
with circulation Tl = 4T/3 and Tr = 2T/3 located at zl(0) = — bo/3 and z#(0) = 26 /3, see Figure [l] 
In both cases, the total circulation is 2r and the separation distance between the two vortices is bo- The 
inviscid (y — 0) evolution of these two cases is well understood. Since the induced velocities from the point 
vortices on each other are always perpendicular to the connecting line, the distance between the vortices 
remains constant bo for all time, and the vortices are co-rotating with a constant rotation rate around the 
origin 

all of which is well known and discussed in [19]. The rotation angle can be integrated from Q. The 
position of the left vortex is zl = — (bo/2)e %e and the right vortex zr = (bo/2)e %e for the symmetric case, 
whereas zl = —(bo/3)e %e and zr = (2bo/3)e ze for the asymmetric case. In general, the inviscid evolution 
of N point vortices is known to exhibit regular, including fixed and relative equilibria, as well as chaotic 
dynamics depending on the number of vortices, their strengths and initial positions, see [TT| \TE\ fl9| [20] . 

Our goal in this paper is to describe dynamical evolution of these co-rotating vortices for non-zero 
viscosity (y ^ 0). Instead of numerically solving the Navier-Stokes equation, we use a simple, analytically- 
tractable, low-dimensional model. For a comparison between the model and Navier-Stokes equation with 
a different initial configuration, the reader is referred to [2 . The model assumes that the vorticity of each 
initial point vortex spreads axisymmetrically as an isolated Lamb-Oseen vortex, thus modeling the diffusion 
term uAuj in 0, while its center moves according to the local velocity induced by the presence of the other 
diffusing vortices, thus accounting for the convection term u- Vcj in ([!]). Stated differently, the model assumes 
that the vorticity field is a superposition of multiple Lamb-Oseen vortices 
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where r = vt. The associated velocity field is computed by substituting ([5| into (J2|, 
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And the velocity of the center of the /3 th vortex is given by subtracting the effect of itself from (J6|, 



The system of equations in ([5|, (J6|, and ^ is referred to as the multi- Gaussian model, or the core-growth 
model, discussed as the basis for a numerical scheme, for example, in [3 j. It is useful to notice in ^ that the 
entire numerator can be thought of as the vortex strength T a (z a , r) which is both time-dependent, and 
depends on the position of each vortex with respect to all of the other positions. Thus it is not surprising 
that the transient dynamics associated with the model depends on the initial configuration chosen. We will 
augment this basic model by coupling it with an appropriately chosen passively advected field, discussed 
further in Section 2. 

1.2 Some relevant literature on vortex merger 

The merging of two co-rotating vortices is an essential phenomenon in transitional and turbulent flows 
and as a result, has emerged as a "canonical" problem in the fluid dynamics literature. It seems to be 
the "simplest" problem whose evolution involves the interplay between two key terms in the equations of 
motion for incompressible flows. First is the convective term given by u • Vcj in the vorticity form of the 
Navier-Stokes equations, second is the viscous term given by uAu. Of course both operate simultaneously 
in evolving the initial configuration, through the transients, to the final Lamb-Oseen state. The convective 
term is usually thought of as being responsible for aspects of merging such as filamentation and stretching 
("inviscid" aspects), while the viscous term is of course responsible for the ultimate diffusive merging. 

Some studies have focused primarily on the inviscid aspects of merging, such as Dritschel & Waugh [6] 
and Trieling et al. j7] who studied both symmetric and asymmetric vortex pairs, and identified 5 flow regimes 
of vortex interaction, namely, elastic interaction, partial straining-out, complete straining-out, partial merger 
and complete merger. Trieling et al. verified their analytical results with experiments, using dye, i.e. passive 
particles, visualization in a viscous flow field. However, in a viscous context, as pointed by Gallay & 
Wayne [1 , the vorticity field will always eventually approach a single Gaussian, hence, the vortex pairs will 
always result in a complete merge. 

Others have focused on vortex interactions in viscous flow. In the case of symmetric co-rotating vortex 
pairs, Melander et al. [8 examined the flow structure in a co-rotating frame, and identified the inner cores, 
the exchange band and the outer recirculation regions. Cerretelli & Williamson [9 and then Meunier et 
al. [10] described the merging process in four stages and they also identified the filaments as the primary 
reason for the merger. However, Velasco Fuentes [11] argued that filaments are not the primary reason for 
the merger, instead they are merely the result of the merger. In the case of asymmetric co-rotating pairs, 
there are fewer studies. Brandt & Nomura [12] classified regimes of interaction based on the time scales, and 
associated the merging with strain rate of velocity field. 

In reality, both the convective and the diffusive terms operate simultaneously. Untangling the "order" 
in which their presence is felt is not a well defined problem, and is analogous to answering a "chicken vs. 
egg" question. Among the many complications in elucidating these issues separately is the fact that they 
depend on the Reynolds number associated with the system. This is most clearly brought out in the paper of 
Cerretelli & Williamson [9] who analyzed the flowfield in a time-dependent co-rotating reference frame and 
were able to convincingly identify four distinct stages of the merging process (their Figure 4) at Reynolds 
number 530: (i) the first diffusive stage; (ii) the convective stage; (iii) the second diffusive stage; and (iv) the 
merged diffusive stage. The first diffusive stage is primarily responsible for spreading the initially confined 
concentrated vorticity profile to an outer "recirculation" region of the flow whose velocity field redistributes 
the vorticity to form two pairs of weaker vortex dipole structures (i.e. a vortex quadrupole) which by its 
nature pushes the original stronger vortices together, when the second and merged diffusive stages act to 
finish job. Our low-dimensional model described in this paper captures well the first diffusive stage, the 
second diffusive stage, and the merged diffusive stage, but cannot generate quadrupole vorticity structures, 
hence does not capture the essential feature of the convective merging stage. By this, we do not mean to 
imply that the model ignores all convective (inviscid) aspects of the transient dynamics, only those that are 
able to generate the quadrupole moment of vorticity. This is the essential deficiency in using the core growth 



(7) 



3 




0.5 1 

(a) vs. r 

Figure 2: (a) Rotation rate and (b) rotation angle for the co-rotating vortex pair. Parameters values are: r = l,6o = 
I, Re = 1000. 

model to study two-vortex merging, but also points to a low-dimensional way to correct the deficiency which 
we discuss at the end. 

Determination of a merging criterion for the vortex pairs has been the focus of many studies both 
in inviscid and viscous setups, mostly in the symmetric case [13j [14j [15j [16], and fewer in the asymmetric 
case [6j[7J[T2]. We study a symmetric case with identical circulations, T, initially assigned to each vortex, 
and an asymmetric case with unequal circulations 4T/3 and 2T/3 initially assigned to each vortex. Both 
cases have the same total circulation of 2r and are shown in Figure [I] Both cases will approach the same 
asymptotic Lamb-Oseen solution in long time since their total circulation is the same. However, the transient 
dynamics is rather distinct and interesting, which will be shown in the later sections. For the symmetric case, 
the onset of merging is usually determined by the ratio between the vortex core size a and the separation 
between the vortex centers 6, and the critical ratio (a/b) cr has been found to be between 0.23 and 0.31 in 
the mentioned studies. For the asymmetric case, whether a universal critical ratio can be found still remains 
an open question. 

The organization of this paper is as follows: we implement the multi- Gaussian model in the context of 
the co-rotating pairs in Section [2] We describe the evolution of vorticity and absolute velocity fields for the 



symmetric and asymmetric cases in Section |2.1| We then show the relative velocity field in Section |2.2[ and 
focus on the evolution and bifurcation of the instantaneous fixed points and instantaneous separatrices. Since 
weaker vorticity tends to be "passively" advected by stronger more concentrated vorticity in the full equations 
of motion (i.e. "one-way" coupling of the vorticity field), a key diagnostic tool for us is the tracking of passive 
tracers that are seeded judiciously in the flow field. The evolution of these passive patches are presented 
and discussed in Section |2.3| Finally, in Section [3] we make connections between the results obtained from 
the multi- Gaussian model and the previous studies and we point to ways in which the "low-dimensional" 
core-growth model used in this paper can be augmented with quadropole terms to increase its fidelity to the 
full Navier-Stokes system. 

2 Viscous evolution of symmetric and asymmetric mergers 

We examine the dynamics of the examples given in Figure [l] using the multi- Gaussian model ([5|,(|6| and 
We non-dimensionalize the problem by length scale L = bo, and time scale T = 27r 2 6o/r which is the inviscid 
turnover time. The Reynolds number in viscous vortex dynamics is defined as Re = Y jv. 



2.1 Vorticity and absolute velocity fields 

Symmetric case For the symmetric co-rotating pair, the vorticity field in ^ takes the form: 
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Figure 3: Symmetric vortex pair: (a) Vorticity along the connecting line. The two vorticity peaks eventually decay and merge 
to a single peak at the origin, (b) Distance between the vorticity peaks and origin as function of time. The peaks merge at 
time r = 0.125. (c) Vorticity evolution at the origin. The maximum vorticity occurs at r = 0.0625. 
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Figure 4: Symmetric vortex pair: Streamlines of the absolute velocity field in inertial frame at r = 0, 0.25, 0.05 and 0.075. 
Each plot shows a [—1, 1] x [—1, 1] window in (x, y) plane. Parameters are r = 1, bo = 1, Re = 1000. 



Due to symmetry, similar to the inviscid case, the vortex centers are also in relative equilibrium, rotating 
around the origin. However, the rotation rate is no longer a constant, instead it is given by 
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Integrating 6 in time, one has the orientation angle 6 as a function of time 
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in which the Euler function is denned as Ei(x) = f^ oo exp(x)/x dx in the sense of principle value. Figure [2] 

shows 6 and 9 as functions of the scaled time r. Initially, 6(0) is the same value as the inviscid case Q, then 
gradually decays to zero as r — > 00. The positions of vortex centers are given by 
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The vorticity along the line that connects the vortex centers is plotted at various instants in Figure |3(a)| 
The distance between the symmetric vorticity peaks and origin is denoted r p . At r p , o ne ha s, duo/dr = 
and d 2 uj/dr 2 < 0, r being the distance from the origin. The result is plotted in Figure 3(b) The turning 
point is obtained when the second derivative is equal to 0, there fore r = 0.125. The vorticity at the origin 
is given by ujo = (T/2ttt) exp(— 6q/16t), and is plotted in Figure 3(c) The maximum vorticity at the origin 
occurs when r = 0.0625, which is obtained by solving r such that duo/dr = 0. 

The absolute velocity field is given by substituting (11) into (Jo]) , which gives 
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Figure 5: Asymmetric vortex pair: (a) Vorticity along the connecting line. One vorticity peak eventually decay and merge 
with the valley, (b) Distance between the vorticity peaks (and valley) and origin as function of time. One peak and the valley 
merge at time r « 0.0725. (c) Vorticity evolution at the origin. The maximum vorticity occurs at r w 0.0304. 




(a) r = (b) r = 0.008 (c) r b w 0.0146 (d) r = 0.05 

Fi gure 6: Asymmetric vortex pair: Streamlines of absolute velocity field in inertial frame at r = 0, 0.008, 0.0146 and 0.05. 
Each plot shows a [—1, 1] x [—1, 1] window in (x, y) plane. Parameters are r = 1, bo = 1, Re = 1000. 



The streamlines associated with the velocity field are depicted in Figure [4] Two main features can be 
observed: (i) the flow field eventually evolves into a single Gaussian vortex as predicted; and (ii) while 
the structure is evolving, it is also rotating unsteadily with rotation rate given b y (|9| ). The instantaneous 



stagnation points in the flow field can be obtained by setting the right hand side of (12) to be zero, and their 
elliptic or hyperbolic characters can be determined by analyzing the eigenvalue problem associated with the 
linearized system around the instantaneous stagnation points (see [2] for a detailed analysis of the collinear 
state). Initially, two elliptic points coincide with the two vortices and a hyperbolic point at the origin. As 
r > 0, the elliptic points begin to merge towards the origin, then at a finite bifurcation time 75 « 0.050, they 
collapse at the origin, transforming the origin into an elliptic point. The topology of streamlines changes 
accordingly. Among the streamlines, the ones that pass through the instantaneous hyperbolic point are called 
instantaneous separatrices, i.e. the thick lines in Figure [4j with arrows representing the directions of the flow. 
Initially, the separatrix is in oo-shape centered at the origin, and the flow field is divided into three regions. 
As time evolves, the separatrix shrinks towards the origin while the elliptic points converge. Eventually, at 
Tfc, the separatrix collapses at the origin, and the flow field is characterized by only one type of streamlines. 
Note that, contrary to a steady state flow field, a time-dependent separatrix does not constitute barriers 
to fluid motion and fluid particles typically move across this separatrix. Interestingly, the bifurcation time 
of the velocity field (75 ~ 0.05) does not coincide with the time that vorticity peaks collapse at the origin 
(r = 0.125). 



Asymmetric case For the asymmetric co-rotating pair, the vorticity field in ([5| takes the form: 
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The rotation rate and angle are the same with the symmetric case given in ([9| and (10). Indeed, the rotation 
rate only depends on the total circulation of the two vortices and their initial distance, not on the strength 
ratio of the two vortices. The positions of vortex centers are given by 

ZL = -\b e^, z R =^b e ie . (14) 



The vorticity along the connecting line for the asymmetric case is plotted at various instants in Figure [5(a/ 



The distances between the vorticity peaks and origin are denoted tl and tr for the respective (initial) left 
and right vortices. The distance between the vorticity valley and origin is denoted as ry. The result is 
plotted in Figure |5(b)[ Similar to the symmetric case, the vorticity field eventually becomes single peaked. 
However, instead of the two peaks and one valley all collapse to form the single peak as in the symmetric 
case, the right peak and valley collapse when r « 0.0725 at vr = ry = bo/2, and annihilate each other, 
and the left peak remains to be the single peak in the field. The evolution of vorticity at origin is ujq = 



(r/67rr)[2exp(-6g/36r)+exp(-6§/9r)], it is plotted in Figure 5(c) The maximum vorticity at origin occurs 
when r « 0.0304. 

The streamlines of asymmetric case in inertial frame are depicted in Figure [6] at various instants. 
Initially, the topology of separatrix is the same oo-shape as that of the symmetric case, though the left and 
right rings are in different sizes. However, at bifurcation time T5 ~ 0.0146, the hyperbolic point collapses 
with one of the elliptic points and they vanish, only one elliptic point remains. The position of the remaining 
elliptic point approaches the origin as time goes to infinity. Note that, similar to the symmetric case, the 
bifurcation time of the velocity field does not coincide with the time that right peak and valley merge in the 
vorticity field. 



2.2 Relative velocity field 

It is convenient to observe the dynamics in a frame co-rotating with the vortex pair with the time-dependent 
rotation rate 0. Let £ = (£, rj) denote position of a point in this rotating frame. The transformation from 
the rotating to the inertial frame is given by 



z = #£, R 
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sin£ 
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(15) 



Symmetric case. For the symmetric case, the positions of vortex centers in this rotating frame remain 
at £l = ( — fro/2,0) and = (6q/2,0) for all time. As time evolves, the vorticity field, initially concentrated 
at £l and begins to spread spatially, inducing a velocity field similar to that of a Rankine vortex with 
time-dependent core R c . The reader is reminded that the velocity field of a Rankine vortex is such that for 
the region inside the core, the velocity field corresponds to a rigid rotation, while for the region outside the 
core, the velocity field decays proportionally to the inverse of the distance. In Figures [7| we superimpose the 




Figure 7: The velocity field induced by the co-rotating vortex pair becomes analogous to that of a Rankine vortex. The 
component of velocity v v along the £ axis is depicted at various instants (solid) . The velocity of a Rankine vortex with vorticity 
2T and time-dependent core are superimposed (dashed). The velocity field is similar to a rigid rotation close to the origin and 
an inverse decay at larger distance from the origin. 
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(a)r = (b) r = r* « 0.068 (c) r = 0.1 (d) r = r 2 * « 0.118 




(e) r = 0.12 (f) t = Tg ~ 0.127 (g) r = 0.13 (h) r = r 4 * « 0.132 

Figure 8: Evolution of separatrices associated with the relative velocity field of the non-symmetric co- rotating pair in rotating 
frame. The arrows on separatrices show velocity direction. Hyperbolic points are represented by cross-sections of separatrices. 
Elliptic points are small circles. 



velocity field of a Rankine vortex on the vertical component of the velocity field induced by the viscously 
evolving co-rotating pair along the horizontal axis at three different instances. The core size R c is obtained 
by equating T/ttR^ = 9, and it increases in time. As the distance from the origin increases, the velocity 
field of the co-rotating pair structure decays analogously to the inverse decay with vorticity 2T. Since we 
are interested in the dynamics close to the origin, motivated by this analogy with the Rankine vortex, we 
examine the time evolution of the relative velocity field obtained by subtracting a rigid body rotation from 
the absolute velocity field expressed in the rotating frame, namely 

£ = R T z - R T RR T z = R T z - 6^ , (16) 

where £ x = (—77, (). Note that the second term corresponds to the rigid body rotation of the rate 0. Hence 
the relative velocity is the difference between the rotated absolute velocity field and a rigid body rotation, i.e. 
the difference between the solid line and the dashed line close to the origin in Figure [7j The instantaneous 
separatrices and stagnation points (elliptic points represented by small circles and hyperbolic points by the 
intersections of the separatrices) associated with (Jl6| are depicted in Figure [8] Initially, one finds a total of 
7 stagnation points: a hyperbolic point at the origin (0,0), a pair of elliptic points on the horizontal axis 
(±6o/2, 0), a pair of hyperbolic points on the horizontal axis denoted by (±£*, 0) (initially £*(0) = boy/5/2), 
and a pair of elliptic points on the vertical axis (0, ±60 \/3/2). The elliptic points on the vertical axis remain 
stationary and of elliptic character for all time. The origin and the pair (±&o/2,0) remain fixed points 
for all time but their characters change when the pair (±£*,0) passes them as it approaches the origin. 
Consequently, the topology of instantaneous separatrices changes in time, as depicted in Figure [8] The 
topological bifurcation times r| ~ 0.118 and r| « 0.132 correspond to when (±£*,0) collide with (±&o/2,0) 
and the origin, respectively, while rj* « 0.068 and r| « 0.127 correspond to the collapse of separatrices 
in this process. After rj, the topology of the separatrices remains the same. The sequence of topological 



bifurcations is summarized in Figure 12 The topology of the separatrices at the bifurcation values are 
depicted in boxes. 

Asymmetric case. Following a similar procedure, we examine the bifurcation sequence of the relative 
velocity field of the asymmetric case. Initially, 7 stagnation points exist: 3 hyperbolic points and 4 elliptic 
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(b) r = 0.03 (c) r = t* « 0.0510 (d) r = 0.05105 




(e) r = r 2 * » 0.0511 (f) r = 0.0515 (g) r = r 3 * » 0.0520 (h) r = r 4 * « 0.0524 
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(1) 



oo 



Figure 9: Evolution of fixed points and separatrices associated with the relative velocity field for the symmetric co-rotating 
pair in rotating frame. The arrows on separatrices show velocity directions. Hyperbolic points are represented by cross-sections 
of separatrices. Elliptic points are small circles. 



points. Notice since the flow field is not symmetric about the vertical axis (but symmetric about the 
horizontal axis), the pair of elliptic points is not located on either axes (the other 5 stagnation points are on 
the horizontal axis). As time evolves, the 1 st , 2 nd and 4 th bifurcations are due to collapsing of separatrices; 
while the 3 rd and 5 th bifurcations are due to the colliding of instantaneous stagnation points; these 5 



bifurcations occur within a very short period of time (r* 



0.0510, r 2 * 



0.0511, T 3 * « 0.0520, r 4 * « 0.0524 



and Tg ~ 0.0525). The topology of separatrices remains the same with the second to last state after the 
5 th bifurcation time. The last bifurcation occurs only when r — >• oo, though the meaning of instantaneous 
stagnation points and separatrices diminishes now since the velocity field is zero everywhere. The topological 



sequence of bifurcations is summarized in Figure 14 



2.3 Passive tracer evolution 

As we mentioned in our previous work [2], passive tracer evolution in the mult i- Gaussian model reveals useful 
information about the vorticity evolution in the Navier-Stokes equations. This is because weaker vorticity 
tends to be passively advected by stronger more concentrated vorticity in the full equations of motion (i.e. 
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(g)r 4 *^ 0.132 



(h) r = 0.16 



(i) r = 0.2 



Fi gure 10 1 Passive tracer evolution for the symmetric case. 



one-way coupling of the vorticity field). Interestingly, the evolution of passive tracers by the fluid velocity 
field is tightly linked to the streamline evolution of the relative velocity field shown in Section 2.2 This 
point is discussed further in Section [3] Since at time r = the vorticity is concentrated at two points, 
we let the vortex cores diffuse for a short amount of time To, then seed the flow with passive tracers of six 
different colors as shown in Figure [To] to distinguish tracers within the left (red) and right (blue) cores of 
radii ao = v^to, left (yellow) and right (teal) rings of inner radii ao and outer radii 2ao, and left (purple) 
and right (green) rings of inner radii 2ao and outer radii 3ao- Since ao is the core, i.e. the standard deviation 
of a Gaussian function, the 3ao range accounts for more than 99% of the total initial vorticity distribution, 
known as the 3-sigma rule. We let the passive tracers be advected by the fluid velocity field given in 
The snapshots of tracer evolution are rotated into the rotating frame for easier comparison. 



Symmetric case. We let To = 0.003, hence ao = ^/^tq ~ 0.110. The initial state is shown in Figure 10 'a). 
As time evolves, the tracers first rotate around the two elliptic points which coincide with the vortex centers. 
Before r = 0.0235, they are confined within the inner separatrix (the region bounded by the inner separatrix 
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(d)r = 0.062 (e)r = 0.1 (f) r = 0.2 

Fi gure 111 Passive tracer evolution for the asymmetric case. 



is referred to as the inner cores hereafter), see Figure 101)). After r = 0.0235, some of the outer tracers 



(purple and green) start to "leak" into the band between inner separatrix and middle separatrix (referred 
to as the exchange band), and they mix while circulating in the band around the inner cores. As the middle 
separatrix continues to shrink, some of the outer tracers will consequently leak outside the exchange band 
when they encounter the hyperbolic points (along the unstable manifolds associated with the hyperbolic 
points). At the first bifurcation time Tf ~ 0.068, the exchange band disappears, the inner tracers still 
remain in the inner cores, some of the outer tracers are in the inner cores while the rest of the outer tracers 
now erode in the regions bounded by the outer separatrix and the middle separatrix (referred to as the outer 
recirculation area) which start to form filaments, see Figure 10 'd). Since the separatrices are time dependent 
Eulerian streamlines, the passive tracers will not faithfully follow the separatrices. Therefore when the tracers 
travel along the outer separatrix and recirculate close to the opposite hyperbolic point, they will "escape" 
to the region outside the outer separatrix (referred to as the outside area), see Figure [lO^e). In the mean 
time, some of the middle (yellow and teal) and later inner (red and blue) tracers encounter the hyperbolic 
points and enter the outer recirculation area, they mix with the outer tracers while traveling close to the 
outer separatrix. At the second bifurcation time r| ~ 0.118, the moving hyperbolic points collide with the 
elliptic points, and most of the inner tracers still remain close to the now hyperbolic points (±&o/2, 0), see 
Figure 101). As the separatrices continue to evolve, the tracers continue to mix in the outside area, but 
closer to the origin the six colored tracers are separated by clear boundaries. After the last bifurcation time 
r| ~ 0.132, the topology of the separatrices remain the same, and the shape of the streamlines will not vary 
much (although the magnitude of the velocity will continue to decay). Hence, after r| ~ 0.132, the tracers 
confined within the inner separatrix (which include all six colors) will mostly remain inside, and rest of the 
tracers will circulate around the outer separatrix in the outside area. As shown in Figure [lOfh) and (i), due 
to the difference in the magnitude of velocities inside the inner separatrix, the boundaries between the inside 
six colored tracers will start to stretch, causing the tracers to mix in layers. 



Asymmetric case. We let To = 0.002, hence ao = y/^Fo ~ 0.089. The initial state is depicted in Fig- 
ure 11 ^a). As time evolves, the inner separatrix shrinks and the green tracers leak into the exchange band 
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Fi gure 12: Homotopic equivalences of the separatrices for symmetric co-rotating pair. Bifurcation states are depicted in 
boxes. 




Fi glire 13: Timeline of important events for the symmetric case. 



and circulate around the inner cores, see Figure [TTJb). As the inner separatrix continue to shrink, all of the 
green tracers will enter the exchange band, followed by all the teal tracers and some of the blue tracers, while 
the left tracers still remain inside the inner core and rotate around the left elliptic point. This is because 
in this case, the right ring of the inner separatrix shrinks, but the left ring actually expands. As the right 
hyperbolic point continue to approach the right elliptic point, some of the tracers leak along the unstable 
manifold. As mentioned before, the 5 bifurcation times occur within a short period of time, hence the po- 



sitions of the passive tracers will not change much during the bifurcation times. As shown in Figure 11 ^c), 
the exchange band disappears, all the left tracers remain in the inner core, while some of the right tracers 
circulate and the others leak to the outer recirculation area or the outer area. The majority of blue tracers 
remain close to the now hyperbolic point. In Figure [Tl^d) and (e), the blue tracers escape along the unstable 
manifolds, and the right tracers erode into the band between the outer and inner separatrices. As the left 
hyperbolic point continues to approach the left elliptic point, the left ring of the inner separatrix continues 
to shrink, and the purple tracers start to leak outside the left ring, see Figure [lljf). Since the left ring will 
shrink to zero only when r — » oo, it will take infinite amount of time for the red tracers to completely leak 
along the unstable manifold associated with the left hyperbolic point. 



3 Discussions and conclusions 

There are three key features in the analysis of the core-growth model we highlight here: 

1. The streamline description of a relative velocity field (velocity field minus a rigid body rotation) and 
their topological bifurcations analyzed in an appropriately chosen time-dependent rotating frame; 

2. The passive trace field description in this same rotating frame; 

3. The identification of the passive field as a "one-way coupled" model of vorticity evolution giving insights 
into the dynamics of a weak vorticity field (the outer region of a decaying Gaussian core) advected by 
the more strongly concentrated vorticity field emerging from the concentrated initial conditions. 

The description offered by the model reveals much more detailed and delicate structure then could 
easily be pinpointed by careful experiments addressing the same questions, such as [9 . Their experiment, 
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Figure 15: Timeline of important events for the asymmetric case. 



however, clearly reveals the importance of the "two-way coupling" of the vorticity and velocity fields in the 
convective stage, as the vorticity is squeezed out of the inner zones shown in our Figure [lOj with respect to the 
generation of a quadrupole structure which is the key "inviscid" mechanism responsible for pushing the two 
vortex centers towards each other. This feature cannot be captured in our "one-way coupled" description, 
although the general monotonically decreasing trend of the distance between vortex centers as shown in our 
Figure [3(b)| is at least consistent with the general direction. 

The rotating frame and topology of streamline patterns is essential to the understanding of viscous 
vorticity evolution. In the multi-Gaussian model, the angular velocity of the structure (which will eventually 
decrease to zero due to the presence of viscosity) is a natural choice of rotation rate, and the analogy between 
the velocity field of the vortex pair and that of a Rankine vortex justifies the physical meaning of subtracting 
a rigid body rotation from the (rotated) absolute velocity field. Indeed, although the passive tracer evolution 
is a result from integrating the absolute velocity field, it is obvious that the relative velocity fields (Figure 8] 
and [9| are much more relevant to the tracer evolution than the absolute velocity fields (Figure [4] and [6). 
The reason was essentially noted by Velasco Fuentes [11] — the Lagrangian dynamics of the passive tracer 
trajectories do not always follow the time varying Eulerian streamlines. In steady flows, both geometries 
are identical. The slower the Eulerian velocity field is varying, the less difference between the two [21]. In 
this case, the absolute velocity field is a fast varying Eulerian field, which is both rotating and evolving. 
The relative velocity field is a much slower varying Eulerian field since it is a result of subtracting the fast 
varying part (rigid body rotation) off the absolute field. 

By using the multi- Gaussian model, we are able to analyze the bifurcation sequence of the relative 
field separatrices, thereby obtaining the topological bifurcation times. The results can be partially verified 
in comparison to the experimental results obtained by previous studies of the symmetric case. The best 
example is Figure 13 in Cerretelli & Williamson [9], where they obtained separatrices topologically identical 
to Figure |8ja) and (h), but the intermediate topological states were not explicitly available, and indeed we 
have not proven that they exist for the full Navier-Stokes equations, only for the multi- Gaussian model. 
For the asymmetric case, to our knowledge, no previous research has shown separatrices in the rotating 
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frame other than the initial state. We present the complete sequence in Figure [8] for the symmetric case 
and Figure [9] for the asymmetric case, and one can see that the bifurcations are rich and not at all trivial. 



The timeline of all important events based on the multi-Gaussian model are plotted in Figure 13 for the 



symmetric case and in Figure 15 for the asymmetric case. 

The path is now clear for a corrected low-dimensional core-growth based model of vortex merger. 
One needs to augment the basic description in this paper with the development and growth of an appropri- 
ately chosen quadrupole structure surrounding the co-rotating state, chosen with time dependent vorticity 
strengths which develop as the primary vorticity from the original core is squeezed out of the inner core 
region. It is interesting to note that some of these multi-pole features are systematically used in the work of 
[22] . Implementing these key features in the context of a low-dimensional model of merger will be the focus 
of a separate paper. 
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